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Especially for our German DUG-Members! 

Besonders fur unsere deutschen DUG-Mitglieder! 

Fur das Friihjahr 1993 planen wir ein Treffen der deutschen Mitglieder der DERIVE USER GROUP. 

Ich schlage als Termin den Samstag, den 24. April 1993 vor. Als Ort der Zusa mm e nk unft hat Hr. Ap- 
pel aus Kronach Schweinfurt angeboten. Hr. Appel hat auch bereits freundlicherweise Vorsorge fur 
Nachtigungsmoglichkeit, ein Tagungslokal und fur die PC-Ausstattung getroffen. Herzlichen Dank 
dafur, lieber Herbert. 

Ich konnte mir vorstellen, dass wir uns am Freitag, den 23. April abends in Schweinfurt treffen und 
dass wir dann am Samstag DERIVE-bezogene Themen behandeln. Bitte teilen Sie mir mit der beige- 
legten Information fur den Herausgeber des DNL mit, ob Sie an diesem Treffen teilnehmen wollen 
(mit oder ohne Nachtigung in Schweinfurt). Teilen Sie mir bitte weiters mit, ob Sie bereit waren, einen 
Kurzvortrag (kurze Themenbeschreibung und Dauer) im Rahmen des 1. Deutschen DUG-Treffens zu 
halten. Als Meldeschluss bitte ich den 31. Janner 1993 zu beachten. Sie werden anschlieBend umge- 
hend liber das Zustandekommen und fiber das Programm verstandigt. Das Treffen in Birmingham im 
Friihjahr war eine wirklich gelungene Veranstaltung. Ich lade Sie herzlich ein, sich am diesem Treffen 
in Schweinfurt zu beteiligen. Natiirlich sind auch unsere Freunde, die nicht in Deutschland wohnhaft 
sind herzlich willkommen, wenn sie die Anreise nicht scheuen. Also bitte. denken Sie an den 31. Jan- 
ner 1993. 


DERIVE - BOOK - SHELF 


Hoogewiijs A, DERIVE, Een Steun bij het Wiskunde-Onderwijs? (DERIVE, A Support in Teaching 
Mathematics?) in Wiskunde en Onderwijs 18 (1992) 


AMMU (Arbeitsgruppe Modemer Mathematik Unterricht - Working Group Modem Mathematics 

Teaching). A group of Austrian teachers deals with Information Technology and Teaching 
Mathematics, Computer Algebra Systems (especially DERIVE), Working with Graphing 
Calculators, Didactics, etc. 

1 st issue in Oct. 1992, 82 pages. 

Contact: Dr. Peter Schuller, A-2351 Wr. Neudorf, Reisenbauerring 6/1/35, AUSTRIA. 

New information for the revised version: 

AMMU website: www . ammu . at 


Download all DNL-DERIVE- and Tl-files from 

http : / / www . austromath . ac . at/ dug/ 
http : / / www . bk- teachware . com/ main . asp?session=37 505 9 
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Lieber Derive Anwender! 

Das zweite Bestandsjahr der DUG ist nun ab- 
gelaufen und wenn ich vor einem Jahr stolz von 
205 Mitgliedem berichten konnte, so sind wir 
nun bereits eine Gemeinschaft von iiber 520 
Derivianem auf der ganzen Welt geworden. 

Fast taglich erreichen mich Anmeldungen und 
Anfragen. Die aktuelle "Verteilungsfunktion" 
finden Sie im Kasten. Ein Blick auf die Beitra- 
ge in diesem DNL, sowie eine kurze Vorschau 
auf die nachsten Beitrage zeigt die erfreuliche 
Internationalist unserer Gruppe: Mechanics 
(Finniand), Extended Fourier Series (Brasilien), 
Physics (Kanada), Analytische Geometrie im R 3 
(Osterreich), Iterationsverfahren (Deutschland), 
Ideen fur die Schule (Hawaii). Aus der CSFR 
und aus Ungam ist Mitarbeit angekiindigt wor- 
den. 


Dear Derive User, 

The DUG has now been existing for two 
years and when I could proudly report in the 
Editor's Letter about 205 members one year ago 
so we have grown up to now to a community of 
more than 520 Derivians from all over the 
world. Nearly every day I receive applications 
or/and inquiries. Y ou will find the actual "dis- 
tribution function" in the box. One glance to the 
contributions of this DNL and a short preview 
to next articles show the pleasing interna- 
tionality of our group: Mechanics (Finland), 
Extended Fourier Series (Brasil), Physics (Ca- 
nada), Iterations (Germany), Analytical Geome- 
try (Austria), Ideas for the Classroom (Hawaii). 
Users from the CSFR and from Hungary have 
announced their cooperation. 


Argentina 1 

Germany 160 

Portugal 4 

Australia 5 

Hongkong 1 

Slovenija 1 

Austria 76 

Ireland 2 

Spain 7 

Belgium 3 

Italy 17 

Sweden 23 

Brasil 2 

Luxemburg 2 

Switzerland 23 

Canada 3 

Mexico 1 

United Kingdom 57 

Danmark 2 

Netherlands 18 

USA 81 

Finland 4 

Norway 6 


France 25 

Poland 1 



Bitte beachten Sie die beiliegenden Informa- 
tionen betreffend die Konferenz in Krems 1 993 
und die Bezugsmoglichkeit des Tagungsbandes 
der Spring School 1992 in Krems. 

Sie finden auch ein Papier fur eine Riickmel- 
dung an mich. Bitte niitzen Sie diese Gelegen- 
heit, damit wir Ihnen mit dem DNL und der 
DUG eine noch bessere Serviceleistung bieten 
konnen. Ich ersuche Sie natiirlich wieder, Bei- 
trage fur den DNL zu liefem oder anzuregen. 

Bitte begleichen Sie rechtzeitig den Mit- 
gliedsbeitrag fur 1993. Der Verwaltungsauf- 
wand ist bereits recht groB und ich danke fur Thr 
Verstandnis im Namen meiner Frau, die mir die 
administrative Arbeit ahnimmt. 

Mit den besten GriiBen und Wiinschen zu 
den Festtagen und zum Jahreswechsel 


Please pay attention to the information con- 
cerning the 2 nd Krems Conference on Ma- 
thematics Didactics and the possibility to order 
the Proceedings of the Spring-school 1 992 in 
Krem s . 

Y ou will find another paper which should be 
an information for me. Please take the opportu- 
nity of a feedback so that we are enabled to 
offer more and better services with DNL and 
DUG. I request for deliverance of contributions 
and for suggestions once more. 

Please settle your dues for 1993 in time. The 
administration is becoming enormous and I 
thank for your comprehension in my wife's 
name for she does all the administration work. 

With best wishes and greetings for 
Christmas and for the New Year 1993 
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The Derive-News-Letter is the 
Bulletin of the Derive-User 
Group. It is published at least 
three times a year with a con- 
tents of 30 pages minimum. The 
goals of the D-N-L are to ena- 
ble the exchange of experiences 
made with Derive as well as to 
create a group to discuss the 
possibilities of new methodical 
and didactic manners in tea- 
ching Mathematics. 

Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone/ FAX: 43- (0) 2275/8207 
e-mail: no j o . boehm0pgv . at 


Contributions : 

Please send all contributions 
to the above address . Non- 
Eng, lish speakers are en- 
couraged to write their contri- 
butions in English to reinforce 
the international touch of the 
D-N-L. It must be said, though, 
that non-English articles will 
be warmly welcomed nonethe- 
less. Your contributions will 
be edited but not assessed. By 
submitting articles the author 
gives his consent for reprin- 
ting it in the D-N-L. The more 
contributions you will send to 
the Editor, the more lively and 
richer in contents the Derive- 
News-Letter will be. 


Preview: (Contributions for the next issues) 

Newton's Chaos; 

Extended Fourier Series; 

Physics in the Classroom 

The Mechanics of Erkki Ahonen; 

The Triangle of PASCAL; 


(will be published March 93) 


Impressum: 

Medieninhaber : DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrif t 
Herausgeber: Mag. Josef Bohm 
Herstellung: Selbstverlag 
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Terence Etchells, Bolton, UK 

Tha nk s for the last edition of the DUG I have enjoyed reading it immensely, especially your ar- 
ticle on Riemann at Random. As you know I have done some work in this area and Paul Drijvers poin- 
ted out that my functions (and presumably yours) only worked for monotonic functions. I decided on 
my return to the U.K. that I would try to write a function that would give true Riemann rectangles. 

Y our article has motivated me to write an article on the functions that I wrote, enclosed with this letter. 
I would be most grateful if you would consider publishing it in a future edition of the DNL. I have 
gone some pain to try to explain how the functions work. I sent a copy of the functions to Paul and he 
wrote back commenting that he did not understand how the routine works although he did find it beau- 
tiful! 

DNL: Your article is published in this issue. So are my comments and your answer. I think its 
a fine example for top - down with DERIVE! 

Helmut Dieck., Essen, Germany 

Inzwischen kann ich einigermaBen mit dem Programm umgehen, aber es sind doch noch ein 
paar Fragen offen geblieben. 

1. Mir fiel auf, dass das Programm die Gleichung SQRT(x) = a nicht nachx auflosen kann, es sei 
denn, man setzt fur a eine konrete Zahl (z.B.: a = 3 + 4 i) ein. Warum? 

2. Zu einer Gleichung 4. Grades fand es "Losungen", die keine sind: siehe die beiden Anlagen A1 
und A2. Habe ich da einen Fehler gemacht? Wenn es ein Fehler des Programms ist, welche Feh- 
ler hat es noch? 

Mr Dieck wants to know (1) why the equation SQRT(x) = a cannot be solved for x and he 
(2) created an equation of degree 4. But solving this equation solutions appeared, which are 
no solutions. See the following example: 


(Z - 

4 - 5i ) 

(z - 6 - 7 i ) 

N 

1 

00 

1 

CD 

H> 

2 

66z 

+ 1420z 

- 2315 - 4i 

3 2 

( 6z - 85z + 184z + 275) 


and after a lot of while solving in Exact Mode (about 1 000 seconds calculation time) he finally 
got 20 lines of solutions and after approximating DERIVE gave: 

z=3+4i; z=7+8i; z=3+4i; z=7+8i 

DNL: For your first question see the following hint please: 

1: Jx - a 

2: “Declare Variable x Complex" 

2 

3 : x = a 


The recent versions of DERIVE give more sophisticated answers: 
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#1: 

SOLVEG/x = a, x) = 

(x = IF(a >0, a )) 

#2: 

a :€ Complex 


#3: 

SOLVE C/x = a, x) = 

' ( 7T 

x = IF < - PHASE (a) < 

l L 2 

#4: 

a :e Real (0, ®>) 


#5: 

SOLVE C/x = a, x) = 

2 

(x = a ) 


Derive 6.10 


DERIVE 6.10 has also problems in factorizing the quartic in Exact Mode. Using a numerical 
algorithm DERIVE is able to solve the equation to our satisfaction: 

4 3 2 3 2 

#3: SOLVECz - 20. z - 66. z + 1420-z - 2315 - i..(24-z - 340-z + 736-z + 1100), z) 

#4: false 

4 3 2 3 2 

#5: S0LVE(z - 20-z - 66-z + 1420. z - 2315 - i..(24-z - 340-z + 736-z + 1100), z) 

#6: z = 2.35776 + 4.99075-i. v z = 6.00924 + S. 64223-1. v z = 3.99075 + 3. 35776. i v z = 7.64223 + 7.00924-i. 

4 3 2 3 2 

#7: NS0LVE(z - 20. z - 66-z + 1420-z - 2315 - i..(24-z - 340-z + 736-z + 1100), z) 

#8: z=8 + 9.Lvz = 6+ 7.Lvz = 4 + 5.Lvz = 2 + 3.i. 


Helmut Scheuermann, Hofheim/Taunus, Germany 

Helmut Scheuermann sent an equation which he wanted to solve with DERIVE 2.06. An equation 
which gave no problems using version 1.55 was not to solve with the later version. 


1 : "Version 1 .55" 

x + LN (x) x + LN (x) 

2 : = 

a b 

x x 

3: "soLve " 

4: x = 1 

5 : LN (x) + x = 0 

6: "soLve in approximate Mode: 

7: x = 0.567143 


1 : "Version 2.05" 

x + LN(x) x + LN (x) 

2 : = 

a b 

x x 

3: "soLve" 

1 

4 : LN (x) = 

0 

b a 

5: x (LN(x)+x)-x (LN(x)+x) = 0 


DNL: I wrote the story of this equation from V. 1 .06 up to V. 2.52. Version 2.10 and later is 
able to solve the equation. It seems that some versions have forgotten a part of their earlier 
knowledge. 
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This is the answer of DERIVE 6.10 in exact and approximate mode in 2005: 


x 

x-e =lvx=lva-b=0 

x = IF(a <0, 0) v x = IF(b <0, 0) v x = IF(a > 0, °o) v x = IF(b > 0, ») 
v x = 0.56714B2986 vx=lva-b=0 


Herbert Appel, Kronach, Germany 

Answer to Mr.Neuraths's letter in DNL#7: 

Tha nk your for the hint in DNL#7. When I said "pupils want to test the teachers competence in using 
the program" I meant that the students try to disturb the lesson by changing some menu-options 
(0-P,M-T,W-D etc.). This is especially then annoying when the correct choice is important for the 
solution of a problem on this level. Another example: When I want to teach to solve an equation step 
by step, it is important that the children are not able to use the command SOLVE. Other menu-points 
are not used in secondary schools at all. For example: Calculus, Declare-3DPlot, Complex Numbers 
and so on. Of course it is a welcome effect that the pupils are free in the selection of the way to solve a 
mathematical problem, but the experience (with other programs) shows that the pupils are less confu- 
sed if only the menu-points needed are shown an the display. Y ours faithfully H. Appel 


Gunter Scheu, Pforzheim, Germany 

The first name of Ch. Goldbach was not Charles but Christian. 


See H.Meschkowski, Mathe- 
matiker Lexikon: 


(I found Charles in a textbook 
dealing with number theory. 
Thank you Gunter for giving 
Goldbach his correct Christian 
name. 


GOLDBACH, CHRISTIAN (18. 3. 1690 - 12. 1. 1764) 

Sohn eines Konigsberger Pastors. Studierte Jura und 
hatte seine Freude an der Mathematik. Auf Reisen durch 
Europa traf er mit fiihrenden Mathematikern seiner Zeit 
zusammen (LEIBNIZ; NIKOLAUS BERNOULLI u.a.) und blieb 
mit vielen von ihnen im Briejwechsel, als er Sekretdr der 
Petersburger Akademie wurde. 

In den Protokollen der Akademie finden sich haufig 
Hinweise auf mathematische Mitteilungen GOLDBACHS, 
z.B. iiber Fragen der Reihenlehre, iiber die Riccatische Dif- 
ferentialgleichung. 

Auch der jetzt herausgegebene Briefwechsel mit EULER 
enthalt Mitteilungen, die von dem mathematischen Ver- 
standnis GOLDBACHS und seinem Ideenreichtum zeugen. 
Besonders wichtig ist die bis heute noch nicht bewiesene 
„Goldbachsche Vermutung”: Jede gerade Zahl ist als 
Summe zweier Primzahlen darstellbar. 
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Message 1864 was entered on 5/14/92 at 12:16 PM. (Read 31 times) 

From JASON to PUBLIC about EXPONENTIAL EXPAND 

I am looking for a way to expand exponentials to exponentials instead of 
sines and cosines. Is there a way? I tried changing the branchs but doen ' t 
help. I need this abilty to calculate Fourier Transforms. 

Thanks Jason B. B. 

P. S. I'm still looking for the difference between Derive VER 2.06 and 
2 . 10 

<<< Reply - see message #1865 >>> 


Message 1865 was entered on 5/14/92 at 3:58 PM. (Read 30 times) 

From JERRY GLYNN to JASON about #1864 / EXPONENTIAL EXPAND 

Try MANAGE -EXPONENT I AL-EXPAND-ENTER and then SIMPLIFY. e A (x+y) becomes 
e A x*e A y . 


Message 1888 was entered on 5/19/92 at 7:00 AM. (Read 29 times) 

From JASON to JERRY GLYNN about #1866 / EXPONENTIAL EXPAND 

Jerry, thank you for the reply. I think I found my problem in that when I 
use exp(#i), it automatically simplifies to sines and cosines. If I use the 
j operator instead of #i then I can simplify down with exponentials. I use 
the manage substitute command to replace j with #i before final simplifica- 
tion. This takes care of my problem. Thanks JBB 

Message 2071 was entered on 6/30/92 at 12:12 PM. (Read 21 times) 

From GREG SMITH to PUBLIC about DERIVE/DR-DOS 

Is anyone running Derive (successfully or unsuccessfully) under DR-DOS 6.0? 
I know there were some compatibility problems with earlier versions, and 
wondered if these had been fixed in the latest DR release. A user on one of 
my other BBSs needs information about this. 

<<< Reply - see message #2082 >>> 


Message 2082 was entered on 7/4/92 at 12:06 PM. (Read 16 times) 

From JIM POTTER to GREG SMITH about #2071 / DERIVE/DR-DOS 

I have had no problems using Derive with DR-DOS 6.0. I have the latest up- 
grade that makes it compatible with Windows (Although that is probably not 
relevant to DERIVE.) . I am using SuperStor 2.0 and QEMM instead of the disk 
compression and memory managers that come with DR-DOS. Good Luck. 

Message 2147 was entered on 7/21/92 at 10:48 PM. (Read 17 times) 

From SOFT WAREHOUSE to PUBLIC about EQUATION SOLVING BY BISECTION 

In approximate mode, DERIVE uses efficient numerical techniques to compute 
an approximate real solution of an equation containing only one variable 
(see Section 4.14 in the DERIVE User Manual) . 

For pedagogical reasons, some users have requested a function that would 
show successively more accurate approximate solutions to an equation. If 
u(a) and u (b) have different signs, the function BISECT (u,x, a, b,n) defined 
as follows approximates to n successively more accurate solutions of the 
equation u(x)=0. Note that u (x) should NOT be an equation itself. 
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BISECT_AUX(u , x, a, b, c) := IF(LIM(u, x, a)-LIM(u, x, c) < 0, [a, c] , [c, b]) 

BISECT (u , x, a, b, n) := ITERATES(BISECT_AUX(u , x, ELEMENT(v,l) , ELEMENT (v , 2) , 
(ELEMENT (v , 1) + ELEMENT(v , 2) )/2) , v, [a, b] , n) 

BISECT(x a 3 - 5 • x a 2 + 2, x, 0, 1, 5) 

0 1 

0.5 1 

0.5 0.75 

0.625 0.75 

0.625 0.6875 

0.65625 0.6875 

Message 2168 was entered on 7/27/92 at 10:52 PM. 

From SOFT WAREHOUSE to PUBLIC about ALGEBRAIC SUBSTITUTION 

Some users have had difficulty doing the following exercise presented page 
56 of the DERIVE User Manual: Use the Manage Substitute command to alge- 
braically substitute e for m c A 2 in the expression 

g m c A 2 

SQRT (m A 2 c A 4 + 1) 

The key to successfully performing this exercise is the need to ALGEBRAI- 
CALLY substitute for m c A 2 since syntactic substitution would not replace 
the m A 2 c A 4 term in the expression. The easiest way to do make this alge- 
braic substitution is to substitute e/c A 2 for m in the original expres- 
sion. This substitution is found by solving the equation e = m c A 2 for 
m. 

To actually make the substitution, highlight the ENTIRE expression and then 
issue the Manage Substitute command. When the command asks for a replace- 
ment for m, enter the value e/c A 2. Simplifying the resulting expression 
should give you the expression #2 as shown on page 56. 


2 

SOLVE (e = m • c , m) = 




Message 2238 was entered on 8/19/92 at 1:47 AM. 

From SOFT WAREHOUSE to PUBLIC about LEAST-SQUARES RESIDUALS 

File RESIDUAL. MTH includes a function that computes residuals for a 2D data 
fit computed by the FIT function. The file also includes an example of 
using RESIDUAL. 
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Advice about least-squares fitting: 

It is generally unwise to have more terms than the least of 
a) he square root of the number of data rows, 
b the square root of Option Precision Digits, 
c) he estimated number of significant digits in the data. 

It is wise to determine these terms incrementally as follows: 

1. Plot the data. 

2. Guess a good form for the dominant term and do a corresponding fit. 

3. Plot the residuals. 

4. Quit if they are as small as your estimated random experimental errors. 

5. Reject the most recent term if the residuals are not significantly smal- 
ler than the original ordinates or the previous residuals. 

6. Guess a good form for the dominant term of the residuals and do a revi- 
sed fit of the original data including this additional term. 

7. Go to step 3. 

You will find this complete file on the disk of the year '92. 

Message 2241 was entered on 8/20/92 at 12:47 AM. (Read 11 times) 

From SOFT WAREHOUSE to PUBLIC about CONVOLUTION INTEGRALS 

At a recent trade show we were asked if Derive could do convolution inte- 
grals. We looked up the definition in an electrical engineering text at the 
adjacent booth, and it can be implemented as follows: 

To convolve an expression u(t) with an expression v(t) over the interval 
a<=t<=b, define 

CONVOLUTION (u, v, t, a, b) := INT (u, LIM (v, t, t-x) , x, a, b) 

Usually, v is exp (t) , in which case we can use the more straightforward 
definition 

CONVOLUTION (u, t, a, b) := INT (u, exp (t-x), x, a, b) 

Message 2267 was entered on 8/26/92 at 8:11 AM. (Read 15 times) 

From JERRY GLYNN to PUBLIC about TEACHING GRAPHING 

In Derive we can plot points in 2 dimensions by plotting [a,b] . With some 
of my students who are beginning graphing I suggest that they plot a number 
of points (somewhat at random) . If I show them [2,3] as an example they 
often plot other points with a pair of positive numbers. I may then ask if 
they can plot a point "over here" pointing to the 2nd quadrant or 3rd or 
4th. If they experiment and make it fine... if not I may a suggestion (try a 
negative number) . I also show them how to clean the screen and how to re- 
trieve earlier points. 


Later I ask if they can plot some points in a row (I may also show them how 
to "back up theit viewing helicopter" by using F10) . With a clean screen I 
suggest they plot x and 2x and 3x and leave them to experiment. Can you 
get a graph to slant this way or that way... 

With a clean screen plot a point (maybe it's a planet) and can you plot a 
path through this planet (trial and error all the way) ...how about 2 pla- 
nets ... 3 ... have you tried to plot xx or xx-3? And on and on. 
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Message 2273 was entered on 8/27/92 at 9:18 PM. (Read 12 times) 

From KAREN to JERRY GLYNN about #2267 / TEACHING ALGEBRA I WITH DERIVE 

Jerry - I'm back on line after being gone for too long. School started ye- 
sterday and I used DERIVE today with my ALGEBRA I students. I have several 
discovery lessons written and plan to let the students discover patterns 
rather than memorize rules. I taught at EXETER computer Conference this 
summer using DERIVE and ALGEBRA MENTOR. It went well and the participants 
were excited about using a symbolic manipulator with young high school 
students . 


Message 2282 was entered on 9/1/92 at 1:56 AM. 

From SOFT WAREHOUSE to PUBLIC about CYLINDRICAL COORDINATES 

File CYLINKIN.MTH defines functions that compute the radial, azimuthal and 
axial speed and acceleration of a point p(t) = [r(t), theta (t), z (t) ] mo- 
ving through radian cylindrical coordinates as a function of time t. This 
file can be regarded as additions to the distributed utility file 

DIF_APPS . MTH . 

You will find this complete file on the disk of the year '92. 


Message 2289 was entered on 9/3/92 at 2:24 AM. (Read 26 times) 

From SOFT WAREHOUSE to PUBLIC about CANTILEVER BEAMS 

File CANTILEV.MTH defines functions that compute the shear, bending moment, 
slope and deflection of a cantilever beam having variable load. Young's mo- 
dulus and moment of inertia. The file also contains examples of their use. 

You will find this complete file on the disk of the year '92. 


Message 2293 was entered on 9/4/92 at 12:01 AM. (Read 25 times) 

From SOET WAREHOUSE to PUBLIC about KINEMATICS IN SPHERICAL COORDINATES 

File SPHERKIN.MTH defines functions that compute the radial, longitudinal 
and colatitudinal speed and acceleration of a point r(t) = [radius, - 
longitude, latitude] moving through radian spherical coordinates as a 
function of time t. This file can be regarded as additions to the distribu- 
ted utility file DIF_APPS .MTH . 

You will find this complete file on the disk of the year '92. 

Message 2295 was entered on 9/4/92 at 8:46 AM. (Read 30 times) 

From BENOITF to JERRY GLYNN about OPTIONS/EXECUTE COMMAND 

I have Derive v2.06 and when I try to execute a DOS command it tells me 
"Cannot find COMMAND.COM". My command.com is defined in COMSPEC.I tried to 
put command.com in my default directory (where is DERIVE.EXE) and I have 
the same problem. What is the problem? 

<<< Reply - see message #2296 >>> 

Message 2296 was entered on 9/4/92 at 9:51 AM. (Read 27 times) 

From GREG SMITH to BENOITF about #2295 / OPTIONS/EXECUTE COMMAND 

If you haven't already, try putting a copy of COMMAND.COM in the 
root directory. Programs SHOULD try to read the COMSPEC environment 
variable rather than assuming the name and location of the command 
processor, but many do not. Looking for COMMAND.COM in \ (root) is 
one of the common older strategies . 






plO 
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Message 2311 was entered on 9/9/92 at 1:16 AM. (Read 23 times) 

From SOFT WAREHOUSE to PUBLIC about EULERIAN ANGLE KINEMATICS 

File EULERANG.MTH defines angular velocity in terms of the Eulerian angles: 
precession, nutation and spin. The file also includes an example. 

You will find this complete file on the disk of the year '92. 

Message 2333 was entered on 9/19/92 at 10:22 PM. (Read 16 times) 

From RANDY to PUBLIC about DERIVE & INEQUALITIES 

Does anyone know how (or if you can) plot linear inequalities or systems of 
linear inequalities on Derive. I'm teaching a business linear algebra cour- 
se and I think it would be helpful. 

<<< Reply - see message #2335 >>> 

Message 2335 was entered on 9/20/92 at 7:45 AM. (Read 16 times) 

From JERRY GLYNN to RANDY about #2333 / DERIVE & INEQUALITIES 

I agree and have requested same for quite a while. Dave Stoutemyer has ex- 
pressed concern that the fill algorithm he would use to indicate which side 
of a line was true might "leak" and fill the entire screen if the line had 
a hole in it (like (x A 2-5x+6) / (x-3) ) . Maybe we all need to ask again and I 
guess you have so thanks. 

<<< Reply - see message #2336 >>> 


Message 2336 was entered on 9/20/92 at 3:32 PM. (Read 18 times) 

From DAN to JERRY GLYNN about #2335 / DERIVE & INEQUALITIES 

I agree also! I know that I have expressed this a few years ago. It would 
be nice and practical from a teacher's point of view to be able to solve 
"Absolute Value" in terms both of < > <= and >= both linear and polynomial 
inequalities as well as rational inequalities. I know, memory is a problem, 
but surely, linear and quadratic inequalities with or without absolute va- 
lues would help the Intermediate Algebra student check or verify solutions 
after a study in "Cut-Points" is presented. I would really like to see 
this implemented in "Derive". 

I am convinced, that Derive is a great program for students and teachers of 
Mathematics. Don't leave home without it! 

<<< Reply - see message #2347 >>> 


Message 2347 was entered on 9/24/92 at 8:37 PM. (Read 12 times) 

From KAREN to DAN about #2336 / DERIVE & INEQUALITIES 

I teach Algebra in a high school using DERIVE. Instead of fill which might 
cause leaks why not have DERIVE shade inequalities the way we teach our 
student - Draw diagonal lines. This eliminates the leaking problem and 
still effectively shades the solution region so a beginning student can vi- 
sualize the "answer". 
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Rational Collocation 

A. Marlewski, Pozna'n 

1 . We deal with the rational function F of the form F(x) = N(x)/D(x), where the nominator 

N(x) = a 0 + ci\X + a 2 x 2 + + a„x n and the denominator D(x) = a n+1 + a n+2 x + + a v _ix v ~ n 2 + 

+ x p 11 1 are the polynomials, the degrees of which (respectively not greater than n and equal to 

p-n- 1) are given in advance. We want to determinate their p coefficients a 0 , a u ,a vA such 

that at every given point (xjjy) there holds the equality 

(1) F(^j)=W 

We suppose Xj ^ xy for j ^ k and we call the system of these equalities as the collocation cond- 
itions), the points (xj,Vj) as the collocation knots. 

The task considered by us is known as the (p,n)- th rational collocation problem or the Pade- collo- 
cation problem, because FI. Pade investigated it in details in 1892. Nota bene, the first results con- 
cerning to this question were obtained some decades before, namely in 1821 by A. Cauchy and in 
1846 by C. Jacobi. 

2. The inquired numbers a (l , a\, , a v . i are involved in the system (1). We replace it by the system 

(2) N(xj) =y i • D(xj ),j = 1 .... p, 

the equations of which are the linear. Our new system is equivalent to the system (1) if every col- 
locating knot (xpVj) is accessible. Let's see the following example to understand what does the ac- 
cessibilty mean. 

3. For p = 3, n = 1 and the knots (1,2), (2,2) and (3,3) our system (2) writes 

a 0 + a\ - 2 a 2 = 2 
a 0 + 2 a\ - 2 a 2 = 4 
a 0 + 3ai - 3a 2 = 0 


so it yields the solution a 0 = -6, a\ = 2 and a 2 = -3, in conseqence we get 


F(x) = 


-6 + 2x 
-3 + x 


The graph of this function does not pass through the point (3,3), the resulted function F is equal to 
2 everywhere but not for the argument x = 3, where F is not definded. The point (3,3) is not ac- 
cessible (for the given knots and the function of the given form). In the discussed case system (2) 
is not equivalent to system (1). 


4. In fact, the Pade collocation problem not always has a solution. It is in contrary to the strictly po- 
lynomial collocation problem, wellknown as the Lagrange or Newton interpolating task. The La- 
grange interpolation is nothing else but the Pade collocation having the denominator equal to 1 , 
and it always has the uniquely solution under the natural assumption on the knots (we only need 
xj ^ xy for j *k). 
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5. In Point 7 you can find the unit COL_PADE . MTH and you can use it during your DERIVE sessi- 
on. This unit offers the function COL_PADE having two arguments. The first argument is the 
2-column matrix of the given knots (xj,yj), the number of the rows is the number of these knots 
and, as above, we denote it by p. The second argument is the degree of the denominator, it should 
be the natural number not greater than p. Then DERIVE transforms the expression 
COL_PADE (d, n) into one of the following results: 


a) the text "The 2nd parameter is too big.", ifp < n, 

b) the text "No (unique?) solution.", if the data do not define the solution (at all or uniquely); 
it takes place e.g.,if 


d = 


1 

2 

3 


1 

2 

3 


and n = 1 


there exists no function of the form 


— — passing through the points (1,1), (2,2) and (3,3), 

a, + x 


c) the vector [N(x),D(x)], if there exists the point x, such that D(xj) = = 0; as above, N and D 
stand for the nominator and denominator of the function F we search, and in Point 3 there are 
listed the sample data invoking the situation here described, 

d) the expression N(x)/D(x), otherwise; this is the rational function having the denominator of 
degree < n, the denominator of degree equal to p-n - 1 and the given values y, at the given ar- 
guments x = Xj with j = 1 ...p. 


6. In Point 4 we have noticed that the Pad, collocation with n =p - 1 turns into the classical poly- 
nomial interpolation. The definition of the function COL_PADE lets to call it in this case, too, so 
we can obtain the polynomial N(x) of degree less than p and fulfilling the collocation condition 
N(xj) = y p j = 1 ...p. 

Let obtain it when the data matrix is made of the knots (-1,5/4), (-1/2, 32/25), (1,1/2), (3/2,0), 
(3,1/8), which are visible in Fig.l. Approximating (up to 6 digits) the expression 
COL_PADE (d, 4) gives 

0.0508333x 4 - 0.0308333x 3 - 0.368958x 2 - 0.344166x + 1.19312. 

Thi s expression differs a little from the result produced by the built-in function FIT, which the ar- 
gument is the matrix of six rows, the last five of them are the same of the matrix d and the first 
one (determinating the independent variable and the form of the approximating function) is 

[x, ax A 4+bx A 3+cx A 2+dx+e] . 


(Comment: DERIVE 6 shows no difference between both results.) 
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7. The unit COL_PADE . MTH: 

COL_PADE.MTH by Adam Marlewski (Poznan,Poland), 1992.07.22 
adapted for DERIVE version 6.10 by Josef Bohm, November 2005 


A_(d, n):=VECT0R 


1 im 

.[x ,y]-»[d|jj.l, dij|2] 


B_(d, n) := VECTOR 


DIM(d) - n - 


3.1 


S_(d, n) :: A_(d, n) •B_(d, n) 


n k 

N0M(d, n) I (S_(d, n)) -x 
k=0 k + 1 


k k 

APPEND(VECTOR(x , k, 0, n) .VECTOR (-y -x , k, 0, DIH(d)-n-2)) , j ,1, DIH(d) 


1 

•d , j, 1, DIM(d) 

j,2 


DEN(d, 


n) 


fDIM(d) - n - 2 

I 

L k=0 


(S_(d, 


k"l 

n)) -x 

n + 2 + k v 


DIM(d) - n - 1 

+ x 


C0L_PADE(d, n) := 

If n > DIM(d) 

"The 2nd parameter is too big." 

If DETCA^Cd, n)) = 0 

"No (unique ?) solution." 

If n(LIM(DEN(d, n), x, dj.j_4.1), j_, 1, DIM(d)) = 0 
[N0H(d, n), DEN(d, n)] 

N0M(d, n)/DEN(d , n) 


-1 



1 32 

2 25 


#7 . d := 


1 


1 

2 


3 


0 


3 


1 

S 


C0I PADE(d , 

h- 1 
II 

2.x - 3 


3 

x - 3 


C0I PADE(d , 

2) = 

2 

4 -x - 16.x + 15 


2 

5-x - 11.x + 12 


C0L_PADE(d , 

3) = 

3 2 

20-x - 36-x - 151-x + 240 


207 - 61 -x 


C0L_PADE(d , 

4) = 

4 3 2 

61-x 37 -x 1771-x 

413-x 1909 

1200 1200 4800 

1200 1600 


Some changes were necessary to make Adam's file from 1992 working with DERIVE 6. So 
Adam had two functions called DENOMINATOR and NUMERATOR, which are now imple- 
mented DERIVE-functions (with another functionality). 
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The next screen dumps show the collocation knots together with the Pade approximations. 




There is a pole, which can be seen by zooming out! See finally the quartic! 
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As you can see documented by the following screeshots I was successful in transferring 
Adam Marlewski 's routine on the Tl. I receive the Pade-functions in the PrgmIO-window and 
can plot them together with the scatter diagram given by the knots. Josef 


[» i glbr afca rdotherlpr gia ! ofc lean Upl 


- 1x2 32X25 


- 1X2 32X25 

■ 

1 1X2 

+ m 

1 1X2 


3X2 0 


3X2 0 


[3 1x8 


[3 1X8 


■pade(m,l) Done 

■pade(m,2) Done 

■ pade(m , 3) Done 


pade<m,4)l 


MAIN RAD AUTO 

FUNC 4/30 



7 *T ) 

x 3 - 3 

4 x 2 - 16 x + 15 


5x 2 - 1 1 • x + 12 

-( 20 -x 3 - 36-x 2 - 151 x + 240) 


61 x- 207 


61x 4 37 -x 3 1771 x 2 413 x . 

1909 

1200 1200 4800 1200 

1600 

MAIN RAD AUTO FUNC 5/30 


I am editing matrix m and then calling program pade(d,n) as Adam did in Derive. The rational 
expressions appear on the PrgmIO Screen. The function found is saved as pf and then be 
stored as one of the y-expressions in the Y=-Editor. I stored the four possible functions and 
plotted them together with the given scatter diagram of knots. 

YT^Yfh^Y fs YTfi^Y 
•~itrol |l/Q|Uar|Find...|Hode| 

7n5 

: Local a_,b_,s_, j, k, nun, den, g 
: Local tl,aus,d_,pr 

s<"2nd parameter too big.", "No Solution" 

}+aus 
: 1-^pr 

jdim<d> II l*g 

: If n>g Then 
: l*tl 
: Goto end 

MAIN FiftD AUTO FUNC 

The right screen shows the start of the Pade-programm. 



i» 1 qebrafca rdotherlprqra ! ofc ! ean Up! 1 

■ pade( 

■ pade 

■ pade 

f>,6) 
1 2 
2 2 
.3 3 
1 2 
2 2 
.3 3 

, 1 
,2 

Done 

Done 

Done 

pade <[[1.2- 

[2,2][3.3]],2> 

MAIN RAD AUTO FUNC 10/30 



61 X- 207 

61x 4 37 x 3 1771 x 2 413 x . 1909 

1200 " 1200 _ 4800 ~ 1200 1600 

No Solution 
2nd parameter too big. 

C2x-6 x - 3> 

*£.^ + 3 

2 2 3 

mhin tan aura rune io/jo 


Here you can see how Adam's routine deals with questionable or „no solutions". This is 
again an example for the fact that ist is not very difficult to exchange ideas from one platform 
to the other. I am no expert on this field, I followed Adam's instructions .... 
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True Riemann Rectangles 

Terence Etchells, Bolton, UK 

In the Derive Newsletter number 7, Josef Bohm described a technique for displaying and calculating 
Riemann rectangles of random width. I had independently written functions that would do the same 
thing (except for random widths) at the Spring School in Krems. Paul Drijvers, Utrecht drew my atten- 
tion to the fact that the routines did not give truly upper or lower Riemann rectangles for non mono- 
tonic functions. 

The screen dump below describes the problem. The function sin x between 0 and 71 for seven rectan- 
gles: 



The rectangles are generated by a VECTOR-command shown below in ASCII-format. (Now it is in 
DERIVE Unicode format.) 

rect(a, b,m,n):=[[a,b;a, b+m] , [a , b+m ; a+n , b+m] , [a+n , b+m ; a+n , b] ] 

rectupper(u,a,b,n):=VECTOR(rect(a+j(b-a)/n,0,MAX(LIM(u,x,a+(b-a) j/n) , 

LIM(u ,x , a+(b-a) ( j+l)/n)) , (b-a)/n) , j , 0 , n-1) 

recti ower(u , a , b , n):=VECTOR(rect(a+j (b-a)/n ,0,MIN(LIM(u,x, a+(b-a) j/n) , 

LIM(u ,x, a+(b-a) ( j+l)/n)) , (b-a)/n) , j , 0 , n-1) 

u is the particular function, a and b are the limits of the particular area and n is the number of rectan- 
gles. The procedure creates matrices that will plot the Riemann upper and lower rectangles. The inter- 
vall [a, b\ is cut into n equal strips, each strip is the base of a Riemann rectangle. In the case of the 
upper Riemann rectangles, the co-ordinates of each side of the strip are substituted into the function u, 
and the maximum of these two values is picked as the height of the rectangle. 
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If there is a turning point within a strip, then this procedure will not pick the true maximum (or mini- 
mum) over the strips interval, Hence as in the picture above the rectangle produced is not a true upper 
Riemann rectangle. I left Krems with the task of writing a procedure, in Derive, that would give the 
true upper Riemann rectangles. I doing so I found that Derive 's ability to symbolically manipulate 
expressions makes it a very special and unique programming language. 


If there is a turning point within a strip then the sign of the derivative at each end of the strip is differ- 
ent. If this is the case we need to find the co-ordinate of the turning point and then substitute this value 
into the function u to get the height of the Riemann rectangle. 

The first thing to do is to determine if there is a change of sign of derivative over a strip. If the interval 
is [a,b] and the interval is split into n strips then the x co-ordinate of each end of the j strip is 


(b-a)j 


and a + 


(b -a) (j + 1) 


we need to substitute these values into the derivative of the func- 


tion. This is done via the limit function, i.e. lim — u . By employing the IF, AND and OR func- 

x ^. a +<t^LL dx 

n 

tions we can determine whether there has been a sign change, i.e. 


IF lim — m< 0 AND lim — u> 0 OR lim — u> 0 AND lim — u < 0, then, else 

, (4-a)y dx , (*-<0(7+1) fix (*— <07 fix . (4-o)tf+ 1) fix 

x->a- 1 li,x x— >a-\ li-x x->a-\ u-x x->a-\ xi-x 


If none of the conditions apply else will be invoked i.e. calculate the height by finding the maximum 
of the heights at each end of the strip. As most strips will not have a turning point, it is important to 
have the else to alleviate unnecessary calculation. If we solve w’(x ) = 0 we use the Newton-Raphson 
method. I have chosen an arbitrary 5 iterations the reader may choose to use more or less iterations. 

I also chose the mid point of the strip as the first approximation, i.e. a + — — ^ ^ + ^ hence the solu- 

2/7 


tion of the equation is found thus, ITERATE x- 


dx (b - a)(2j + 1) . 

— ,x,a + - -,5 will give the x co- 

d- 2« 

— , u 


ordinate of the turning point. This now needs to be substituted into «(x) to get height of the rectangle. 
This expression for the heigth will be lim u . 


It is possible that the turning point is not necessarily the maximum height in the strip, so to be safe use 
the maximum (or minimum) of the end points and turning point. 


The then argument will become 


MAX 


u, lim u, lim u, 


( b-a)(j+\ ) 
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The full functions are displayed how they are edited. I find that I prefer to construct functions within 
the algebra window using F3 and F4. Flowever, functions of this size cannot be edited on the author 
line. So I had to employ a line editor and load the file through Transfer Load. Flowever, the debugging 
routines employed when entering an authored expression do not apply. This makes the debugging 
process doubly hard! Any suggestions Soft Warehouse! 


(This paragraph concerning editing and debugging is - fortunately enough - historic, Josef) 

riemupsum(u,a,b,n) :=VECTOR(rect(a + j • (b - a)/n, 0, IF((LIM(DIF(u , x) , x, 
a + j-(b - a)/n, 0) < 0 a LIM(DIF(u, x) , x, a + (j + 1) • (b - a)/n, 0) > 

0) v (LIM(DIF(u, x) , x, a + j • (b - a)/n, 0) > 0 a LIM(DIF(u, x) , x, a + 

(j + 1) • (b - a)/n , 0) < 0), MAX(LIM(u, x, ITERATE(x - DIF(u, x)/DIF(u, x, 
2), x, a + (b - a) • (2 • j + l)/(2-n), 5), 0), LIM(u, x, a + (b - a)-j/n), 

LIM(u, x, a + (b - a) • ( j + l)/n)), MAX(LIM(u, x, a + (b - a)-j/n), LIM(u, 

x, a + (b - a) • ( j + l)/n))) , (b - a)/n) , j, 0, n - 1) 


riemlosum(u,a,b,n) :=VECTOR(rect(a + j • (b - a)/n, 0, IF((LIM(DIF(u , x) , x, 
a + j-(b - a)/n, 0) < 0 a LIM(DIF(u, x) , x, a + (j + 1) • (b - a)/n, 0) > 
0) v (LIM(DIF(u, x) , x, a + j • (b - a)/n, 0) > 0 a LIM(DIF(u, x) , x, a + 
(j + 1) • (b - a)/n , 0) < 0), MIN(LIM(u, x, ITERATE(x - DIF(u, x)/DIF(u, x, 
2), x, a + (b - a) • (2 • j + l)/(2-n), 5), 0), LIM(u, x, a + (b - a)-j/n), 
LIM(u, x, a + (b - a) • ( j + l)/n)), MIN(LIM(u, x, a + (b - a)-j/n), LIM(u, 
x, a + (b - a) • ( j + l)/n))) , (b - a)/n) , j, 0, n - 1) 


To prove that these functions do work below is a Derive screen dump using these function on sin x 
between 0 and n for 7 rectangles. 



It is possible that the reader may be able to find a function where the Newton-Raphson method di- 
verges from the solution of u\x) = 0, hence not giving the true Riemann rectangles. I would be grate- 
ful to receive any such functions or any other correspondence on these routines (in English please) 
from Derive users. 


File RiemTerence . dfw 
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True Riemann Rectangles 
Part 2 

Josef Bohm, Wilrmla 


At first I tried to make the files from above better “readable”. I followed Temece's ideas and used my 
own work with upper sums and lower sums. 1 believe that the file RIEMANN 2 . MTH will explain itself. 

1: [a :=, b :=, F(x) :=] 

d b - a 

2: Fl(t) := lim — F(x), D(n) := 

x->t dx n 

3: [XL(n) := a + (k - l)-D(n), XR(n) := a + k-D(n), XM(n) := a + (k - 0.5)-D(n)] 

4: [ FMAX(n) := MAX(F(XL(n)) , F(XR(n))), FMIN(n) := MIN(F(XL(n)) , F(XR(n)))] 

5: RELEXT(n) := lim F(x) 

x->ITERATE(x - 3(F(x) , x)/3(F(x), x, 2), x, XM(n) , 5) 

FRMAX(n) := 

If Fl(XL(n)) -Fl(XRCn)) < 0 

6: MAX(RELEXT(n) , FMAX(n)) 

FMAX(n) 

FRMIN(n) := 

If Fl(XL(n)) -Fl(XRCn)) < 0 

7: MIN(RELEXT(n) , FMIN(n)) 

FMIN(n) 


8: “Uppersum - Obersumme” 

T XL(n) 0 

! XL(n) FRMAX(n) 

9: OSUM(n) := VECTOR , k, 1, n 

XR(n) FRMAX(n) 

1 XR(n) 0 

10: “Lowersum - Untersumme” 

T XL(n) 0 

i XL(n) FRMIN(n) 

11: USUM(n) := VECTOR , k, 1, n 

XR(n) FRMIN(n) 

1 XR(n) 0 

Y ou have to author the function f ( x ) : = and the intervall a : = and b : = (Can be done in one vector - 
between [..], see DNL#7). Then plot OSUM (n) or USUM (n) - with n being the number of rectangles. 
Don't forget to set “Approximate before Plotting” in the 2D-Graph Window. 






p20 


Josef Bohm: True Riemann Rectangles Part 2 


D-N-L#8 


12: [ F(x) := SIN(8-x) + 1.5, a := 0, b := 4] 
13: [ F(x) , OSUM(IO) ] 

Figure 1 seems to be ok, but ... 

14: [ F(x) , 0SUM(5) ] 

Figure 2 is obviously no upper sum! 



- 0 . 5 • x 

15: F(x) := SIN(6-x) -e +1.5 

16: [ F(x) , USUM(10) ] 

Figure 3 is not so bad! 





- 0 . 8 • x 

17: F(x) := SIN(10-x)-e + 1.5 

18: [ F(x) , USUM(4) ] 

Figure 4 is wrong 

19: [ F(x) , USUM(2) , 0SUM(2)] 

Figure 5 is wrong 
20: b := n 

21: [ F(x) , USUM(10) , OSUM(IO)] 
but Figure 6 and 7 are correct! 



F igure 4 


Figure 5 
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It is obvious that the condition to look for a turning point will fail, if the slopes in the end points of a 
strip containing such a turning point have the same sign or if there are more extremas within one strip. 
I produced some plots and sent them to Terence an .... I got an answer which will form the next chap- 
ter of our RIEMANN-adventure. 

True Riemann Rectangles Part 3 

Terence Etchells, Bolton, UK 


Dear Joseph, 

Thank you for your letter and screen dumps of my Riemann routines running on 
those functions. You are quite right, if the width of the rectangle contains more than one 
turning point then the routine breaks down. This is because only the end points are sampled. 

When one looks at the function the highest (or lowest) point is found by sampling many 
points in between the end points. You indeed set me a task, but 1 have solved the problem (it 
took all day Sunday to write!). Please find enclosed a disk with a new file RIEMANN3.MTH 
(to save you typing), it is 720 K MSDOS format. The function RIEMUPSUM(u,a,b,n,s) will 
find the true upper sum rectangles (OBERSUM), the lowersum is easily produced by editing 
this function. 

It differs from the original in that it samples s points between the end points, if there is a 
change of sign of f (x) between any sample points as the first approximation. The maximum 
of these points is used. As a consequence the computing times are very long! 60+ seconds on 
a 386 machine. 

rect(a . b . m . n) : = [ [ a , b ; a, b+m] , [a, b+m; a+n , b+m] , [a+n , b+m ; a+n , b] ] 

ri emupsum(u , a, b , n , s) : = VECTOR(RECT(a+j*(b-a) / n , 0 , MAX(VECTOR 
(IF(LIM(DIF(u , x) , x , a+(b-a)*j/n+(a+(b-a)*( j+l)/n-(a+(b-a)*j/n))*k/ s)>=0 AND 
LIM(DIF(u , x) , x , a+(b-a)*j/n+(k+l)*(a+(b-a)*( j+l)/n-(a+(b-a)* j/n))/ s)<=0 OR 
LIM(DIF(u , x) , x , a+(b-a)*j/n+(a+(b-a)*( j+l)/n-(a+(b-a)*j/n))*k/ s)<=0 AND 
LIM(DIF(u , x) , x , a+(b-a)*j/n+(a+(b-a)*( j+l)/n-(a+(b-a)*j / n))*(k+l)/s)>=0 , 

LIM(u , x , ITERATE (x-DIF(u , x)/DIF(u , x , 2) , x , (a+(b-a)* j / n+(b-a)*(2*k+l)/ (2*s*n)) , 5)) , 
Max(LIM(u,x,a+(b-a)*j/n) , LIM(u , x , a+(b-a)*(j+l)/n))) ,k,0,s-l)) , (b-a)/n) , j ,0,n-l) 

You will find this file under the name RIEMANN3 . MTH on your disk of the year together with the 
RIEMLOSUM-function. It was easy to edit and to save as a .MTH file with WORD5. 

Having loaded RIEMANN3 in DERIVE I tried to plot the figures from Part 2 once more. See the re- 
sults (I don't list the functions): 
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-0 . 8x 

RIEMUPSUM(SIN(8 • x) + 1.5, 0, 4, 5, 4) RIEMLOSUM(SIN(10x)e +1.5,0,77,4,3) 



-0 . 8x -0 . 8x 

RIEMUPSUM(SIN(10x)e +1.5,0,77,2,5) RIEMLOSUM(SIN(10x)e +1.5,0,77,2,5) 




Congratulations, Terence!! 

In times of DERIVE 5 and DERIVE 6 when DERIVE has become programmable, we can try 
to make such routines better readable and understandable for everybody. I am using an aux- 
iliary function mysol, which looks for possible zeros of the 1 st derivative in the intervals. 


rupsum(u, a, b, n, s := 0.1, d, u2, i, tp, lh, str) := 

Prog 

d := (b - a)/n 
u2 := 3(u, x) 
i := 0 
lh := [] 

str := [] 

Loop 

If i = n 

RETURN [u, str] 

tp := mysol(u2, x, a + i-d, a + (i + l)-d, s) 
h := MAX(VECTOR(u, x, APPEND([a + t d , a + (i + l).d], tp))) 

str := APPEND(str, [[a + i-d, 0; a + 4 ■ d , h; a + (i + 1) ■ d , h; a + (i + l)*d, 0]]) 

i :+ 1 
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- 0.6-x 

rupsum(SIN(10 -x) -e + 2, 0, 10, 6, 0.2) 



- 0.6-x 

rlosum(SIN(10.x).e + 2, 0, 10, 10, 0.2) 



I have to admit that Terence's routine is much faster than my program. I will try to improve it 
and also try to produce shaded rectangles. 
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Esperimenti didattici con un 
prograirana di Computer Algebra 
Schede di Lavoro 

Rosamaria Castelletti, Segrate 

In the last DNL I gave a summary of an article, first published in "L'insegnamento della matematica e 
delle scienze integrate" by Centro Ugo Morin, 31010 Paderno del Grappa {Italy). Here are some of 
Mrs Castelletti's worksheets announced in DNL#7 : 

TAVOLA 3 

Liceo Scientifico Statale "N. Machiavelli" - Pioltello 

a . s . 1989/90 

DERIVE 

Classe prima 


Con l'aiuto della grafica risolviamo alcune disequazioni di vario tipo individuando le soluzioni nel 
piano cartesiano. 

Prima di tutto entriamo in Derive e, al solito, suddividiamo lo schermo in due fmestre: 


a) Consideriamo la disequazione fratta 


x 1 —x-6 


Per risolverla tracciamo il grafico delle due funzioni 


N : y-x 2 - x-6 D: y = x — \ 

(conviene cambiare scala: Scale/x-scale = 3, y-scale = 5) 


Cercate ora di iindividuare dal grafico gli intervalli in cui le funzioni sono entrambe positive (y>0) 

oppure entrambe negative (y<0) perche in questi intervalli la frazione N/D sara 

Quali sono le soluzioni della disequazione? 

Se volete trovare i punti di intersezione di una curva con l'asse x basta usare il 
comando soLve per l'equazione della curva esaminata cioe quella evidenziata dal 
riquadro bianco 


2x - 3x + 1 

b) Ripete i passi precedenti confratta < 0 

(x - 2)(x - 1) 

Se volete cancellare i grafici precedenti, usate Delete. Se volete lasciare l'ulti- 
mo grafico, usate Butlast. 


c) Vediamo ora di risolvere un sistema di disequazioni 


jx 2 -3x-10<0 
jx-3 > 0 


Questa volta dobbiamo visualizzare dove le disequazioni sono entrambe verificate ("l'intersezzione"), 
sempre dal grafico delle curve. 

Soluzioni: 


x - 2 > 0 

d) Proviamo ora con il sistema < x 2 - 9x + 20 > 0 

x 2 - 6x - 7 < 0 


Soluzioni: 







Rosamaria Castelletti: Worksheets - Schede di Lavoro 


D 


-N-L#8 


p25 


e) Risolviamo infine, sempre tracciando i grafici, una disequazione irrazionale: 

2x < 3^4x 2 -13x + 3 . 

Come prime tracciate i grafici delle due funzioni (p e r la radice potete usare ( ) A (l/2) op- 

pure S Q RT( ) ) 

Soluzioni: 


TAVOLA 4 

Liceo Scientifico Statale "N. Machiavelli" - Pioltello 

a . s . 1989/90 

DERIVE 

Classe prima 

1) Entrate in Derive e tracciate il grafico della retta y - -^x + 5 

II punto di intersezione con l'asse x e A( ) 

II punto di intersezione con l'asse y eP( ) 

2) Tracciate il grafico della retta y = -2x + 5, sullo stesso piano della prima retta. 

II punto di intersezione con l'asse ieB( ) 

Il punto di intersezione con l'asse y e 

3) Dimostrare che il triangolo ABP e rettangolo 

(usare l'inverso del teorema di Pitagora e la formula per calcolare la distanza di 
due punti) 

4) Possiamo affermare che "condizione necessaria e sufficiente affmche due rette siano perpendicola- 
ri e che il prodotto dei loro coefficients angolari sia uguale a -1"? 


TAVOLA 5 

Liceo Scientifico Statale "N. Machiavelli" - Pioltello 

a . s . 1989/90 

DERIVE 

Classe prima 


a) y = 3kx + l 

Incominciamo a considerare alcuni casi particolari. Sostituiamo a k il valore 1 e tracciamo il grafico 
della curva. Ripetiamo assegnando a k i valori 2 e 3. 

E' giusta l'affermazione "l'equazione y = 3kx + l rappresenta un fascio proprio di rette"? 

Quale e il centra di tale fascio? 

b) Consideriamo ora l'equazione y = 3k x + 1 (1) 

Se y = 0, si ottiene l'equazione di primo grado (2) 

Le intersezioni con l'asse x dell'equazione (1) sono le soluzioni della (2) ed il loro valore dipender a 
dai valori attribuiti al parametro k. 
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TAVOLA 7 

Liceo Scientifico Statale "N. Machiavelli" - Pioltello 

a . s . 1989/90 

DERIVE 

Classe terza 

OBIETTIVO: scoprire le proprieta di alcune funzione e curve 

LE RAPPRE SENT AZIONI GRAFICHE DEVONO ESSERE TUTTE RICOSTRUITE SUL FOGLIO DI LAVORO. 
INDICARE CON PRECISIONE I PUNTI NOTEVOLI 

Per gli esercizi 1-5 scegliere il valore 5 nell'opzione Scale prima di tracciare i grafici. 


1) Tracciare sullo stesso piano cartesiano le rappresentazioni grafiche di 
y = ^J 9-x 1 e y = -^9-x 2 3 4 5 

OSSERV AZIONI: 

2) Considerare il fascio di rette di equazione 
y = k(x- 5) 

Evidenziare il parametro k. Usare il comando MANAGE SUBSTITUTE. Indicare 3/4 come valore da 
sostituire a k. Ripetere indicando il valore -3/4. Tracciare le due rette cosi ottenute nello stesso piano 
che contiene le curve del punto 1 . 

OSSERV AZIONI: 


3) Cosa succede s eke minore di -3/4 oppure k e maggiore di 3/4? 

E se k e intemo all'intervallo (-3/4, 3/4)? 

Esplorare cosa succede sostituendo a k diversi valori (comando MANAGE SUBSTITUTE). 

Riportare sul foglio di lavoro il grafico ottenuto, scrivendo accanto ad ogni retta il valore di k che l'ha 
generata. 


4) Tracciare la rappresentazione grafica di 
y = \l 3x-x 2 e y = -\l 3x-x 2 

La curva complessiva ottenuta e 

Ripeti il lavoro dei punti 1 -3 con il fascio di rette di equazione 


y = k 


3 " 

2y ' 


OSSERV AZIONI: 


5) Tracciare la rappresentazione grafica di 
y = -2 + V 4-x 2 e y = -2- V 4-x 2 


OSSERV AZIONI: 
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6) Chiudere la fmestra grafica. Considerare la funzione di equazione 


y = — sin(x) + 2 cos(x) 


2 


Tracciare la rappresentazione grafica della funzione (dopo aver riportato a 1 i valori nell'opzione Sca- 
le). 

La funzione e periodica? Se si, quale potrebbe essere il periodo? 

La funzione ha punti di massimo e di minimo? Se si, sapresti determinare con Derive il valore minimo 
e il valore massimo? 

Descrivi il procedimento che hai usato. 


7) Usate Derive per tracciare il grafico della funzione 
y = x sin(x) 

Tracciare sullo stesso piano le rette y = x e y = - x. 

OSSERVAZIONI: 

8) Tracciate ora la rappresentazione grafica di 
v = x 2 sin(x) 

Nell'esercisio 7 abbiamo tracciato le rette y = x e y = -x. In questo caso quali curve conviene traccia- 
re? 

OSSERVAZIONI: 


9) Esplorate infine la seguente funzione: 
y = x + sin(x) 

OSSERVAZIONI: 
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This is a letter too long to be a part of the DERIVE USER FORUM. It is a lucky chance that I had 
some ideas in the same direction when I worked through the messages of the Bulletin Board Sendee. 
This might be an answer and completion to Mr Robin 's letter (Editor). Read now Mr Robin 's letter: 


I was interested to read in DNL#7, in the Bulletin Board Section 
(incidentally how does one access this facility?) the correspondence 
between Andromeda and Boom-Boom concerning the Normal probability 
integral. (DNL#7, page 7) 

Using DERIVE in Exact Mode NORMAL (1) simplifies to: 


NORMAL(l) = 



+ 1 


2 


although it takes 9 seconds to do so. To do inverse problems, the 
equation say NORMAL (x) = 0.95 is solved in Approximate Mode to give 
1.64489, taking 17 seconds. 

I came across the approximation given for Q(x) in the manual for the 
Casio Programmable Graphical calculator fx7000G. Whilst this was 
useful and quick to use (instantaneous), I found that I could not 
dispose the Normal tables unless I also had an inverse for this 
function. This involves solving the equation Q(x) = a for values of 
a less than 0.5. Newton-Raphson method was very suitable, as Q' (x) 
is of course known, but it was necessary to have a good starting va- 
lue, for the method to work efficently. To do this I found a simple 
approximation for Q(x) namely; 

EXP (-0.43 - 0.45 (x + 0.78) A 2 ) 

Now this can easily be inversed to give: 

x = INVQ(a) = V ( (-LN (a) - 0.43) /0 . 45) - 0.78 

which gives us a good starting point for the Newton-Raphson procedu- 
re. The direct program for finding Q(x) is put into a subroutine. 

The fx7000G takes about 2 seconds to find INVQ(0.05) = 1.644854019. 

I am enclosing copies of the calculator programs for finding Q (prog 
7) and INVQ (prog 9) . I am sure any reader proficient with program- 
ming can easily convert for use in any computer or programmable cal- 
culator. I am also enclosing copies of graphs showing the absolute 
and fractional error in using this crude quadratic approximation. 

Yours sincerely 

You can find informations about the Bulletin Board in your DERIVE 
manual, look under Additional Resources. 

I don't know which computer you are using, but when I approximate 
NORMAL (1) it only needs 0.1 seconds and for solving NORMAL (x) = 0.95 
my machine needs 0.8 seconds. I can't copy the programs and graphs 
mentioned above, because of their print quality , and I don 't have 
the special calculator . Maybe Mr. Robin will send a better copy. 

Thank you. It might be that the following routine will suit you (and 
other readers) more. Please read a part of a letter to D . Stoutemyer 
concerning the same problem: 

I've found the approximating algorithm for the inverse distribution function. It might be interesting 

for some users. 
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Approximation of the Normal Distribution function 
Abramowitz/Stegun: Handbook of Mathematical Funtions 
Dover, New York 1965 

Wahrscheinlichkeitsrechnung , Statistik 

Dietmar Herrman 

Vieweg, Braunschweig 1983 

#1: [p := 0.23164, a := 0.31538, b := -0.35655, c := 1.7814] 

#2: [d := -1.8212, e := 1.3302] 

1 

#3: T(x) := 

1 + P- I x | 

2 

1 - x /2 

#4 : F (x) := e 

7(2-tt) 


#5: 4>(x) := 1 - F(x)-(a-T(x) + b-T(x) + c-T(x) + d-T(x) + e-T(x) ) 

DI5T(x) := 

If x > 0 
#6: tji(x) 

1 - 4>(x) 


The inverse Distribution Function: 



#8: [c2 := 0.010328, cl := 0.80285, cO := 2.5155] 

#9: [dl := 1.4327, d2 := 0.18926, d3 := 0.001308] 

2 

#10: Z(w) := c2-R(w) + cl-R(w) + cO 

3 2 

#11: N(w) := d3-R(w) +d2-R(w) +dl-R(w)+l 

Z(w) 

#12: INV_<t<w) := R(w) 

N(w) 

INV_DIST(w) := 

If w > 0.5 
#13 : INV_<t>(w) 

- INV_0(1 - w) 

#14: [DIST(O) , DIST(l) , DIST(2.97), DIST(-1.5), DIST(-1.5) + DIST(1.5)] 

#15: [0.4999, 0.8413, 0.9985, 0.06680, 1] 

#16: [INV_DIST(0.9985) , INV_DIST(0. 0668) ] 

#17: [2.968, -1.500] 


Plotting the absolute error: 
| NORMAL (x) - DIST(x) | 
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Here are some graphics from Mr Setif. in the first functions he demonstrates 
the Polar Coordinates: (Set coordinate system to Polar!) 


picture 1 


VECTOR 

f k "I 

2-C0S(x) + , k, 0, 10 

, VECTOR 

f k V 

- 2-C0S(x) + , k, 0, 10 

. 

L 5 J 


l 5 J. 


picture 2; #5 in black and #6 in red 


VECTOR • 


C0S(3-x) 


2 



k, 12 


( C0S(3-x) k 1 

VECTOR + , k, 0.5. 11.5 

l 2 4 J 




picture 3 


picture 4 


2-C0S(4-x) + O.S 
2-C0S(6-x) + O.S 


c 

VECTOR 


C0S(12-x) 


3 


k, 1. 8 
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picture 5 


2 -COS 


S-x 


+ 0.5 


2 J 
0 <= x <= 4pi 
' 5-x ^ 


COS 


v 


+ 2 


0 <= x <= 8pi 

( 5-x ^ 
2-COS 

l 8 J 

0 <= x <= 16 pi 



Change the coordinate system to Rectangular 




seti f_7 


VECTOR 

f 

POLYG 

r k ^ ^ 

2-C0S(3-x) + , x, 0, 2 -jt, 15 , k, 0, 8 


S. 

l 4 J ) 
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seti’f_8 



r 

r k ) 

1 

VECTOR 

POLYG 

2-C0S(3-x) + — , x, 0, 2-n, 30 

, k, 0, 6 


l 

L 3 J 

J 


r 

r k ) 

1 "1 

VECTOR 

POLYG 

2-C0S(3-x) + — , x, 0, 2-n, 30 

, k, — , 6 


l 

L 3 J 

2 J 


setif_9 


P0LYG(2-C0S(x) +1, x, 0, 2-n, 7) 
P0LYG(2-C0S(x) + 0.2, x, 0, 2-n, 13) 



f 

f B ' x 1 


- 

POLYG 

2-COS 


+ 0.4 , x , 0, S • JT , 6 


k. 

L 4 J 


J 


f 

f 5 ‘ x 1 



POLYG 

2-COS 


, x, 0, 24 • jt , 72 



k. 

l 12 J 

J 





See some more pictures on the next page. 







p34 


Mr Setif's TREASURE BOX 


D-N-L#8 




setif_l 


POLYG 


POLYG 


POLYG 


POLYG 


POLYG 


f COS(12-x) 


+ 2 , x , 0, 2-rr, 20 


13. x 

2 - C0S| |, x, 0, 24-rr, 48 

12 


2 -COS 


29 -x 
30 
6 


, x, 0, 60-rr, 180 


2. 5- COS 


5 -COS 


x 

I , x i 0 , 14 • rr , 

ft-K 


52 


, x, 0, 14-rr, 52 


POLYG 


POLYG 


POLYG 


Mr) 

W-T-) 


6-x 


SIN 


7 

6-x "|2 


x, 0, 14-rr, 78 


, x, 


0, 14-rr, 78 j 


6- COS 


6-x 


7 J 


x, 0, 14-rr, 78 


SIN 


6-x 



